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Resumen 


Este artículo presenta el procedimiento y la aplicación de análisis de confiabilidad 
con diversos métodos en un problema sencillo, como es la capacidad de descarga 
de una alcantarilla. En el análisis se considera que las tres variables de rugosidad, 
diámetro y pendiente poseen distribuciones de probabilidad normales y se aplican 
para efectos de comparación los métodos de integración directa; el de Monte 
Carlo; el de la transformada de Mellin, con la estimación de primer orden de la 
variancia; el método de Rosenblueth y de Harr del punto de estimación; el método 
del primer orden del segundo momento estadístico (MFOSM), y el método 
avanzado del primer orden del segundo momento estadístico (AFOSM), con dos 
versiones: el de Hasofer-Lind y el de Tang. Se adopta como verdadero el resultado 
de la probabilidad de falla obtenida con el método de integración directa (lo cual 
no es posible en análisis complejos, sobre todo cuando las incertidumbres sean 
importantes) y se comparan los resultados. La mayoría de los métodos puede 
aplicarse y es posible obtener estimaciones confiables cuando los métodos se 
aproximan a comportamientos lineales, pero cuando se tiene no linealidad de las 
variables o bien cuando las incertidumbres se incrementan significativamente, la 
precisión de algunos métodos se deteriora rápidamente. Tal es el caso del método 
del primer orden del segundo momento estadístico. Para métodos con muestras de 
gran tamaño en sus variables, el método de Monte Carlo es el de mayor aplicación, 
pero la confiabilidad del método converge cuando se tiene un gran número de 
simulaciones y no se conoce estrictamente el resultado final de la probabilidad de 
falla; otra limitante importante es que el número de variables puede hacer que el 
problema no tenga una solución práctica. Métodos en los que se emplea el punto 
de estimación (Rosenblueth y Harr) pueden ser muy atractivos desde el punto de 
vista computacional, en la medida en que el número de variables se incrementa y 
pueden parecer muy buenos en su empleo, ya que ofrecen resultados parecidos a los 
obtenidos con el método de Monte Carlo y el de integración directa; sin embargo, 
en caso de que las incertidumbres sean importantes, pueden existir diferencias 
significativas. El método del primer orden del segundo momento estadístico es 
aplicable sólo en casos muy sencillos en los que la función de comportamiento está 
claramente definida y existe una linealidad en las variables; empero, en problemas 
complejos pierde precisión rápidamente. El método del primer orden del segundo 
momento estadístico es muy aplicable y puede tomar en cuenta incertidumbres en 
caso de que el analista decida hacer correlaciones de las variables que intervienen 
en el problema y que la mayor parte de las veces se asocian con incertidumbres; ésta 
parece ser una gran ventaja sobre los demás métodos, ya que es posible involucrar 
variables que en muchas ocasiones se ignoran o deprecian por no poder analizarlas. 
El método de Hasofer en problemas simples parece ser bastante apropiado, sin 
embargo, el método de Tang resulta muy atractivo, al analizar el estado límite de 
falla del problema estudiado. 
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Introducción 


El diseño y análisis de la ingeniería de recursos 
hidráulicos trata con la ocurrencia del agua 
en varios sistemas y sus efectos en diversos 
aspectos, tales como el ambiente, la ecología y 
lo relativo a la sociedad. 

Debido a la naturaleza extremadamente 
compleja de los procesos físicos, químicos, 
biológicos y sociales, se han realizado enormes 
esfuerzos por diversos investigadores para 
tratar de tener un mejor entendimiento de 
los mismos. Un producto benéfico de estos 
estudios de investigación ha sido el desarrollo 
de modelos que describen la interrelación 
e interacción que tienen cada una de las 
componentes de dichos procesos, lo cual se 
pretende en este artículo. El término “modelo” 
se refiere a cualquier elemento estructural y 
no estructural de transformación que produce 
alguna clase de efecto en la salida de los 
modelos mencionados. 

En la ingeniería de los recursos hidráulicos, 
la mayoría de los modelos son estructurales 
y toman la forma de modelos matemáticos, 
ecuaciones, cuadros, gráficas o programas de 
computadora. El modelo es una herramienta 
útil para los ingenieros, que les permite 
predecir el comportamiento del sistema 
ante varios escenarios posibles en los que la 
eficiencia y efectividad del modelo puede 
formular dicho comportamiento. Aun así, se 
debe hacer un mayor esfuerzo para mejorar 
nuestro entendimiento de los procesos que 
intervienen en los sistemas hidráulicos y acotar 
las incertidumbres que se tienen acerca de las 
variables que aparecen en dichos modelos. 

En general, la incertidumbre debida a 
procesos aleatorios estocásticos no puede ser 
eliminada, pero puede reducirse a través de un 
manejo cuidadoso en la recolección y el manejo 
de los datos que intervienen en los diversos 
procesos analizados. 

Las incertidumbres que intervienen en sus 
diversos procesos se dividen básicamente en 
cuatro categorías: hidrológicas, hidráulicas, 
estructurales y económicas. Más específica- 


mente, en los análisis y diseños se tienen como 
origen de las incertidumbres las naturales, las 
inherentes a los propios modelos, las asociadas 
con los datos y las operacionales. 

Las naturales están asociadas con su propia 
naturaleza aleatoria, tales como la ocurrencia 
de precipitaciones y avenidas; los eventos 
hidrológicos tienen variación en el tiempo y el 
espacio, y su naturaleza no puede ser predicha 
con exactitud. Debido a lo anterior, el modelo 
no es más que una abstracción de la realidad, 
que generalmente involucra simplificaciones 
e idealizaciones. La incertidumbre de los 
modelos refleja la discrepancia que se tiene entre 
la realidad y el modelo. Las incertidumbres 
paramétricas son el resultado del manejo de 
datos que miden los fenómenos; también 
pueden tener el efecto de la recolección 
propia de los mismos y, en ocasiones, de 
cambios operacionales que afectan la correcta 
interpretación de los datos. 

Las incertidumbres en los datos pueden 
deberse a (1) (2) 
inconsistencia y no homogeneidad de los 


errores de medición, 
mismos, (3) errores en el manejo y la transcrip- 
ción y (4) una inadecuada representación de 
las muestras analizadas. 

En general, las incertidumbres operacio- 
nales incluyen asociadas 


aquellas con 


la construcción, manufactura, deterioro, 
deficiencias en el mantenimiento y, en muchas 
ocasiones, el error humano. La magnitud 
de este tipo de incertidumbres depende 
fundamentalmente de la correcta operación 
de los sistemas y del control de calidad 
adecuado en la construcción y fabricación de 
los equipos con los que se operan los sistemas. 
El deterioro progresivo debido a la falta de 
un mantenimiento adecuado puede resultar 
en cambios importantes en la resistencia y en 
ocasiones en una disminución de la capacidad 
estructural, lo cual debe merecer un especial 
cuidado en el proceso de análisis de diversos 
sistemas de ingeniería en el futuro. 

El objetivo de este artículo es mostrar 
el estado del arte de cómo evaluar las 


incertidumbres y la confiabilidad en los 
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proyectos hidráulicos de ingeniería, con la 
aplicación a un problema sencillo, como es la 
capacidad de una alcantarilla. 


Técnicas para análisis de incertidumbres 


Hay varias técnicas con diversos métodos 
analíticos que permiten la derivación exacta de 
las distribuciones de probabilidad acumulada 
(DPA) o de los momentos estadísticos de 
variables aleatorias como una función de varias 
variables aleatorias. 

Sin embargo, el éxito de su implementación 
depende fuertemente de la relación funcional, 
formas de la DPA estudiada y herramientas 
matemáticas aplicadas. 

Los métodos analíticos son herramientas 
poderosas en problemas sencillos, aunque su 
utilidad se limita en los problemas complejos 
de la vida práctica (real). Las situaciones que 
existen en los que las técnicas analíticas pueden 
aplicarse para obtener incertidumbres de 
modelos de salida con derivaciones analíticas 
es virtualmente imposible; sin embargo, es 
práctico entonces encontrar una solución 
numérica aproximada. 

La evaluación de incertidumbres y sus 
parámetros básicos, media y desviación 
estándar (dos primeros momentos) es funda- 
mental, debido a la necesidad de conocer 
con exactitud estos valores, que permitirán 
determinar la confiabilidad del sistema de 
ingeniería analizado. 


Métodos para evaluar riesgos 
a) Método del periodo de retorno 


El método del periodo de retorno es un 
procedimiento tradicionalmente empleado en 
sistemas de ingeniería asociados con eventos 
como los hidrológicos y sismológicos, donde 
los fenómenos que producen las cargas oO 
demandas son variables aleatorias producidas 
por fenómenos naturales a lo largo del tiempo y 
entonces se asocia el concepto de ocurrencia o 
recurrencia. De esta manera se fija un periodo 


de retorno de diseño relacionado con eventos 
tales como avenidas, lluvias o sismos que afectan 
las presas, vertedores, puentes, alcantarillas, 
etcétera, obteniéndose un gasto pico o un sismo 
(aceleración del terreno) de diseño; aunque 
puede estimarse la probabilidad de falla y la 
de seguridad, y por tanto la confiabilidad del 
sistema, los periodos de retorno se asocian con 
eventos históricos presentados y se extrapolan 
para estimar el evento de diseño, por lo que su 
aplicación es muy limitada; sin embargo sigue 
siendo muy usada en la práctica profesional. 
En cuanto al método del periodo de 
retorno, Wood (1977) estimó los riesgos al 
evaluarlos con un método de transformación 
(1981) 
consideraron factores de incertidumbre para 


integral; Duckstein y  Borgardi 
integrar la unión de la función de densidad 
de probabilidad de “resistencia y carga”, 
para evaluar la probabilidad de la falla de 
un sistema de presas; Cornell (1969) aplicó 
el método del valor medio del primer orden 
del segundo momento estadístico (MFOSM) 
en la ingeniería de diseño; Tung y Mays 
(1977) aplicaron el MFOSM para estimar las 
confiabilidades estáticas y dependientes de 
tiempo en sistemas de alcantarillado; Wood 
(1977) expandió 
miento del punto de falla utilizando series 


la función de comporta- 


de Taylor, usando el método avanzado 
de primer orden del segundo momento 
estadístico (AFOSM); McKay et al. (1979) 
establecieron el muestreo del hipercubo 
latino (LHS) para mejorar la convergencia de 
la función de comportamiento del método de 
Monte Carlo. Marengo (2006), por su parte, 
empleó el AFOSM para estimar el riesgo de 
falla en obras de desvío, en los que comparó 
la probabilidad de falla obtenida con este 
método y el del periodo de retorno. 

Este periodo de retorno se selecciona por 
guías y recomendaciones; por ejemplo, la 
avenida de diseño para la obra de excedencias 
es de 10000 años para puentes de 200 0 250 años; 
para la obra de desvío en cortinas flexibles 
(tierra y enrocamiento) de 50 a 100 años, y para 
cortinas de concreto (arco y gravedad) de 10 o 
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20 años; para el análisis sísmico de la cortina 
de una presa se considera un sismo de diseño 
con periodo de retorno de 1 000 a 10 000 años; 
se estima como sismo más frecuente el asociado 
con 200 años de periodo de retorno. 

El periodo de retorno se define como el 
tiempo promedio en que una magnitud de la 
resistencia Y será igualada o excedida (Chow, 
1953; Ang et al., 1991). Por lo anterior, si T, se 
expresa en años (para T_ > 1), la probabilidad 
de que un evento Z iguale o exceda a Y en cada 
año está dado por: 

jl 


P(Z>Y)= =— 


TV) (1) 


En fenómenos naturales, Z se supone como 
una variable continua, y si el riesgo de falla se 
define como la probabilidad de ocurrencia de 
que Z sea mayor que Y en cada año, entonces 
la probabilidad de que un sistema no falle es: 


1 (2) 
TO 


P(Z<Y)=1- 


La probabilidad de que no se presente una 
falla en alguno de los n años de la vida útil del 
proyecto será: 


pan=| -=m) 


1,0) dd 


Cabe mencionar que al desarrollar las 
ecuaciones (1), (2) y (3) se han hecho dos supo- 
siciones: a) la ocurrencia de los fenómenos de 
la variable aleatoria Z es independiente en la 
vida útil de n años, y b) el sistema natural es 
invariante en el tiempo. 

Sin embargo, cuando un sistema hidráulico 
y, en general una estructura, queda expuesto 
a variaciones temporales, entonces la 
probabilidad asociada con el periodo de 
retorno, como se expresó en las ecuaciones 
anteriores, tampoco puede usarse como una 
medida del riesgo de estructuras sujeta a este 


tipo de acciones (Borgman, 1963). 


Aun cuando el riesgo evaluado con la 
ecuación (3) considera una buena parte de 
los riesgos por eventos naturales, existen 
incertidumbres asociadas con las variables 
que integran las expresiones de carga y 
resistencia de los sistemas analizados que no 
se toman en cuenta, por lo que el riesgo total 
de un sistema complejo en el que intervienen 
fenómenos naturales no debe evaluarse con 
otro método. 


b) Método de integración directa 


El riesgo se evalúa por medio de una integra- 
ción directa, analítica o numérica de las 
funciones de densidad de la carga y resistencia. 
En este caso, las funciones de distribución 
mencionadas deben estar definidas. Si dichas 
las 
variables que representan, entonces el método 
es exacto. 


situaciones describen correctamente 


En 1980, Tang presentó un procedimiento 
que permite incorporar las incertidumbres 
al modelo de probabilidad en la evaluación 
del riesgo de falla en presas; para ello 
utilizó una integración directa para la 
valoración del riesgo hidrológico. En 1977, 
Wood calculó la sobreelevación y el riesgo 
estructural en forma analítica, con funciones 
de densidad de probabilidad supuestas 
para las avenidas y para los modos de falla, 
donde no se consideraron las incertidumbres 
hidrológicas, hidráulicas y de modelación 
paramétrica. También en 1977, Tung y Mays 
definieron los riesgos para alcantarillas y 
bordos, al estimar primero los parámetros 
estadísticos de carga y resistencia a partir 
de las incertidumbres de los parámetros con 
fórmulas de aproximación de primer orden 
y asignando entonces funciones de distribu- 
ción a la carga y resistencia. 

El riesgo evaluado es sumamente sensible a 
la función de distribución asignada, por lo que 
una suposición impropia o una aproximación 
mal hecha de dichas acciones pueden deme- 
ritar la exactitud ganada al hacer la integración 
directa. 
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Duckstein y Bogardi, en 1981, también 
estudiaron el problema de bordos asociados 
con varios tipos de falla, como sobreelevación, 
tubificación, deslizamiento de taludes y erosión 
por viento. El riesgo se estimó con la integra- 
ción directa de la unión de las funciones de 
densidad de las variables de carga y resistencia; 
sin embargo, la selección de la resistencia fue 
un tanto ambigua y las incertidumbres no se 
representaron claramente. 

La mayor desventaja del método de 
integración directa es la gran dificultad 
que se tiene para derivar correctamente las 
funciones de distribución de probabilidad 
de las variables de suministro y demanda, 
especialmente en sistemas complejos, como 
las presas de tierra y enrocamiento. Asi- 
mismo, una vez que las funciones de densidad 
se establecen, existe una gran dificultad 
para integrarlas, incluso con la ayuda de 
computadoras. De esta manera, el método 
de integración directa es bueno solamente 
para sistemas simples o cuando se requiere 
gran exactitud en la evaluación del riesgo y 
se conocen perfectamente las funciones de 
densidad de probabilidad de las variables 
que intervienen en el problema. 


c) Método de simulación de Monte Carlo 


Es un proceso que utiliza en cada simulación 
un conjunto particular de valores de variables 
aleatorias generadas artificialmente de 
acuerdo con la distribución de probabilidad 
que se está analizando. 

Con cálculos relativamente sencillos, pero 
repetitivos, se puede encontrar un conjunto 
de valores con los que se estima la función de 
distribución de probabilidad del margen de 
seguridad F,,(0). El riesgo de falla se evalúa 
como la relación entre el número de valores 
positivos o negativos de F,,(0) con respecto al 
total de elementos generados. 

En 1981, Duckstein y Bogardi estimaron 
la probabilidad de falla en presas de jales 
(minas). En 1972, Haan evaluó con este 


método las probabilidades de error en 


problemas hidrológicos, en función del 
número de observaciones usadas para la 
determinación de parámetros de los modelos 
estocásticos. 

En 1975, Matalas et al. aplicaron el método 
para estimar los parámetros estadísticos de 
media, desviación estándar y coeficientes de 
asimetría a varias distribuciones de secuencias 
de flujo. 

Chow, en 1953, utilizó el método para 
generar secuencias de datos para el estudio 
del comportamiento de sistemas hidrológicos. 
En 1977, Wen usó los resultados obtenidos con 
la simulación de Monte Carlo para verificar la 
derivación de la estadística de combinaciones 
de cargas extremas. 

Ang et al. (1991) manejaron el método 
de simulación de Monte Carlo con el 
llamado muestreo de importancia en casos 
multidimensionales; Leira (1991) lo aplicó 
con distribuciones normales multivariadas, 
e Ibrahim y Rahman (1991) lo usaron para 
revisar la confiabilidad de sistemas dinámicos 
con incertidumbre. 

Según Hoshiya et al. (1991), dicho muestreo 
de importancia parte de considerar que la 
estimación de la probabilidad de falla de 
diversos sistemas estructurales se hace en un 
rango tan pequeño (10* o 107), que la función de 
densidad de probabilidad es postulada como 
una función desconocida del tipo exponencial 
F(Z, a,, a, ...) en la vecindad de Z = 0, donde 
A,, A, ... SON Valores constantes a identificarse 
(generalmente por correlación). 

En la medida que el muestreo de Z, se 
inicia, se obtiene inmediatamente un grupo 
de datos experimentales para la función de 
distribución de Z usando Z, y se hace un ajuste 
de los datos con los que se consigue F(Z, a,, a, 
...). Aunque este ajuste es inestable en un 
principio, en la medida que se procede en el 
muestreo, el proceso es esencialmente secuen- 
cial y se puede llevar a cabo por medio del pro- 
ceso de filtrado de Kalman, el cual se describe 
claramente en la referencia de Hoshiya et al. 

Cuando a, converge a valores en los que el 
coeficiente de autocorrelación es cercano a uno, 
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la probabilidad de falla del sistema se estima 
como: 
P,=P(Z<0)=F(0, di, Apr...) (4) 

Cabe subrayar que con este método es 
posible obtener una convergencia eficiente 
debido a que el procedimiento es una 
“aproximación lineal” de P,, en lugar de una 
estimación puntual, como tradicionalmente se 
hace. 

De hecho, puede decirse que la técnica 
de Monte Carlo es tal vez la única solución 
técnica a problemas que no pueden resolverse 
analíticamente debido al comportamiento 
no lineal o complejo de las relaciones que 
intervienen en los sistemas analizados. Sin 
embargo, este método tiene las siguientes 
desventajas: 


1. El riesgo estimado al usar la técnica no es 
único, depende del tamaño de la muestra y 
del número de simulaciones. Los momentos 
estadísticos reales de la unión de funciones 
de distribuciones de probabilidad no son 
del todo ciertas. 

2. El costo y tiempo de computación que se 
consume con esta técnica se incrementa 
sustancialmente en la medida que el nivel 
de precisión y el número de variables se 
incrementen, en general, se recomienda que 
si se puede aplicar un modelo analítico, éste 
debe preferirse a la simulación de Monte 
Carlo. 


Confiabilidad 


La confiabilidad de un sistema de ingeniería 
“es más realista cuando se mide en términos 
de la probabilidad” (Tang, 1984). El objetivo 
de un análisis de confiabilidad (Marengo, 
2006) es asegurar durante la vida útil de un 
proyecto o en el periodo de tiempo en el que es 
evaluado que X > Y, donde Xes la capacidad 
de resistencia del sistema y Y es la capacidad 
de demanda (o carga); colocados tanto la 
demanda como la resistencia en un entorno 


natural, ambos están sujetos a variabilidad 
por diversas condiciones de carga externas; 
es obvio que también la resistencia puede 
cambiar en el tiempo. 

Tradicionalmente en un problema de 
(Marengo, 2006), la 
confiabilidad se expresa en función del factor 


demanda-resistencia 


de seguridad FS = X/Y o como el margen de 
seguridad MS = X-Y, mientras que las varia- 
bles FS, MS, Xy Y se consideran determinísticas. 
Si la resistencia y la demanda tienen naturaleza 
estocástica, también FS y MS serán variables 
aleatorias estocásticas. 

Cuando un análisis se lleva a cabo con 
variables estocásticas (Marengo, 2006), los 
resultados usualmente se expresan en términos 
del índice de confiabilidad f, al que se le 
calcula la probabilidad de ocurrencia p, que en 
un sistema de ingeniería hidráulica se define 
como la probabilidad de seguridad (no falla) en 
el que la resistencia del sistema excede la carga; 
esto es: 

p.=P(L<R) (5) 

La probabilidad de falla, P, es el 
complemento de la confiabilidad o riesgo, que 
puede expresarse como: 


R=P,=P(L>R)=1-P, (6) 

Debe recordarse que el riesgo en un siste- 
ma de ingeniería puede estimarse (USACOE, 
1976) como: 


1 
R=1-(1--)" 
( —? 


F 


(7) 


Donde R = riesgo, T.= periodo de retorno 
y n = el periodo de análisis del evento (que 
usualmente es la vida útil), ya que es posible 
asociar la probabilidad de ocurrencia con el 
periodo de retorno (ecuación (3)); entonces, 
cuando se emplea este criterio, el riesgo se 
expresa como: 


R=1-F, (8) 
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Hay dos enfoques básicos probabilísticos 
para evaluar la confiabilidad de un sistema 
de ingeniería hidráulica. El enfoque más 
directo es un análisis estadístico de los datos 
de los últimos registros de falla para sistemas 
similares. El otro enfoque es a través del 
análisis de confiabilidad, que considera y 
combina la contribución de cada factor que 
potencialmente influye en la falla. El primero 
es un enfoque de sistema “lumped”, en el que 
no hay necesidad de ningún conocimiento 
sobre el comportamiento de la historicidad de 
la estructura ni su carga ni su resistencia. Por 
ejemplo, los datos de falla por desbordamiento 
en presas revelan que el mayor valor 
promedio ocurre para presas entre 15 y 30 m 
de altura, con un 50.82% de los casos (21 de 
61 casos de fallas totales) (Marengo, 1996) y la 
probabilidad de falla para presas construidas 
en nuestros días es alrededor de 10* por presa 
por año (Marengo, 1996). En muchos casos, 
este enfoque directo no es práctico, porque 
(a) el tamaño de la muestra es demasiado 
pequeño para ser estadísticamente confiable, 
especialmente para la baja probabilidad de 
consecuencia de eventos; (b) la muestra puede 
no ser representativo de la estructura o de la 
población, o de grupos específicos de la misma; 
y (c) las condiciones físicas de la presa pueden 
ser no estacionarias, es decir, que varían con 
respecto al tiempo. El riesgo promedio de falla 
de presa antes mencionado sólo considera los 
casos de falla totales y no distingue presas 
viejas de nuevas presas, en las que las prácticas 
de seguridad reducirán seguramente el riesgo 
de falla; tampoco distingue el tipo de falla que 
se está analizando. 

Hay dos pasos importantes en el análisis 
de confiabilidad: (a) identificar y analizar 
la incertidumbre de contribución de cada 
factor y (b) combinar las incertidumbres de 
los factores estocásticos para determinar 
la confiabilidad general de la estructura. El 
segundo paso, a su vez, podrá proceder de 
dos maneras: (1) directamente, combinando 
las incertidumbres de todos los factores 
o (2) las 


separadamente, al combinar 


incertidumbres de los factores pertenecientes 
a diferentes componentes o subsistemas 
para evaluar primero la confiabilidad de 
subsistema respectivos y, a continuación, 
confiabilidades de 
diferentes componentes o subsistemas a 
ceder la capacidad de confiabilidad global 
de la estructura. La primera forma se 


combinando las los 


aplica a estructuras muy simples, mientras 
que la segunda es más conveniente para 
sistemas complicados. Por ejemplo, para 
evaluar la confiabilidad de una presa, las 
confiabilidades asociadas con la hidrología, 
hidráulica, geotecnia, estructuras y otras 
disciplinas, podrían evaluarse por separado 
primero y, a continuación, se combinan para 
encontrar la confiabilidad general. O bien, 
de los componentes, podían ser evaluadas 
en primer lugar, según los diferentes modos 
de falla y, a continuación, combinarlos. 
Vrijling (1993) proporciona un ejemplo real 
de la determinación y combinación de las 
confiabilidades de las componentes en el 
diseño de la barrera que surge de tormentas 
en la zona oriental de Holanda. 

En esta 


sección se describen varios 


métodos analíticos que permitirían una 
derivación exacta de la DPA y/o momentos 
estadísticos de una variable aleatoria como 
una función de varias variables aleatorias. 
El éxito de la implementación de estos 
procedimientos depende en gran medida de 
la relación funcional, las formas de las DPAs 
involucradas y la habilidad matemática del 
analista. Los métodos analíticos son potentes 
herramientas para problemas que no son 


demasiado complejos. 
a) Técnica analítica: “Transformada de Mellin” 
Cuando las variables aleatorias en una función 


W = g(X) son independientes y no negativas, y 
g(X) tiene una forma multiplicativa como: 


N 
W=3(X)=2.[[ Xx” (9) 


iZ 
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la transformada de Mellin es especialmente 
atractiva para realizar análisis de incertidumbre 
(Tung, 1990). Para una función con DPA f(x), 
donde x es positiva, se define como: 


M,(s)=ML£(0)] 


a (10) 
=Í, cdi ()dx=E(x%), x>0 

Donde M,(s) es la transformada de Mellin 
de la función f(x) (Springer 1979). Por lo 
tanto, la transformada de Mellin proporciona 
una forma alternativa para encontrar los 
momentos de cualquier orden para variables 
aleatorias no negativas. 

Al igual que la propiedad de convolución 
de las transformadas Exponencial y de Fourier, 
la transformada de Mellin de la convolución 
de las DPA asociadas con múltiples variables 
aleatorias independientes en una forma de 
producto es simplemente igual al producto 
de las transformadas de Mellin de las DPA 
individuales. Además de la propiedad de 
convolución, la transformada de Mellin tiene 
varias propiedades operacionales útiles, como 
se muestra en el anexo 1. 


b) Técnica aproximación: método de estimación 
de primer orden de la variancia (FOVE) 


Este método es llamado también método 
de la propagación de la variancia. Estima 
incertidumbres en un modelo que se basa 
en propiedades estadísticas de modelos de 
variables aleatorias. La idea básica del método 
es aproximar un modelo que involucra varia- 
bles aleatorias por medio de la expansión de 
series de Taylor. 

Considerando que una función hidráulica o 
hidrológica W está relacionada con N aleatorias 
X¡, Xy «Xy COMO: 

W=gX)=2(X, X, --. Xy) (11) 

Donde X = (X, X,, ..., X,)' es un vector 
columna N dimensional de variables aleatorias 
y t = representa la matriz transpuesta. 


La expansión de la serie de Taylor de la 
función g(X) con respecto a las medias de las 
variables aleatorias de X = Hu en el espacio 
puede expresarse como: 


| (Xu) 


Wes 


(12) 


Na (a%(x) 
+ A) (UOH) +8 

Donde H, es la media de las variables X, , 
es la media de las variables X, y e representa los 
términos de orden superior. 

Las derivadas parciales de primer orden 
son los llamados coeficiente de sensibilidad y 
cada uno representa la razón de cambio en el 
modelo W con respecto a la unidad de cambio 
de cada variable en hu. 
de orden 
superior (e en la ecuación (12)), la esperanza 


Despreciando los términos 


del modelo W puede expresarse como: 


N 


1 3 (9%00 
E(W)= 8(1) El Cov[X;, X¡](13) 


Y la variancia de W= g(X) puede expresarse 
como: 


Var [W]= Al al E 40% -1)] 


il ja 


N N N 2 
E) y dex) [0 2(X) (14) 
2444 oX; le 9X; 0X; q 
El (XX, 0% 110) ] 
Cuando las variables aleatorias están 


correlacionadas, la estimación de la variancia 
W, usando la aproximación de segundo orden, 
requeriría del conocimiento del producto 
de momentos cruzados entre las variables 
correlacionadas. En la práctica, esta infor- 
mación raramente está disponible; cuando las 
variables son independientes, las ecuaciones 
(13) y (14) pueden plantearse como: 
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y: 
N a 2 
var (25) 2 
j (16) 
1% /|0e00| |0%(x) 
pele] 
i=l p Xi p 


En esta ecuación (16), la variancia de W, a 
partir de la aproximación de segundo orden, 
bajo la condición de que todas las variables son 
estadísticamente independientes, requeriría del 
conocimiento del tercer momento. 

En las aplicaciones prácticas, donde los 
momentos de orden superior y los productos 
cruzados de los momentos no están disponibles 
fácilmente, la aproximación de primer orden 
frecuentemente se adopta. 

Al truncar los términos de segundo orden y 
mayores en la serie de Taylor, la aproximación 
de primer orden de Wen X= hu es: 


E(W) =3(1) =W (17) 


Var [W]= A a 03 (X) 
ia ja Xx; 0X; m (18) 


Cov(X,, X;) 


y =s! C(X)s 


En la cual s = V_ W(u) es un vector columna 
N-dimensional de coeficientes de sensibilidad 
evaluados en el punto de falla; C(X) es la matriz 
de variancia-covariancia del vector aleatorio 
X. Cuando todas las variables aleatorias son 


independientes, la variancia del modelo W de 
salida puede aproximarse como: 


Var(W)= ys 0 =D (19) 


i=l 


En la que «G representa la desviación 


A pon : 2 2 2 
estándar y D = diag(o,?, 0,, ..., 0,?) es una 


matriz diagonal que involucra las variancias 
de las variables aletorias. De la ecuación (19), 

Siof 
Var(W) 
incertidumbre en el modelo de salida que 


la relación indica la proporción de la 


contribuye con la incertidumbre de la variable 
aleatoria X.. 


c) Técnica aproximada: método de Rosenblueth 
del método de estimación probabilística (EP) 


El método de estimación de aproximación 
probabilística de Rosenblueth (EP) es una 
técnica computacional directa que permite 
de 
puede usar para estimar estadísticamente 


hacer un análisis incertidumbre. Se 
momentos de cualquier orden, en un modelo 
que 
correlacionadas o no. 


involucra varias variables aleatorias 


Originalmente, el método de Rosenblueth 
fue desarrollado para usarse con variables 
sintéticas (Rosenblueth, 1975). Posteriormente 
se extendió para el tratamiento de variables no 
sintéticas (Rosenblueth, 1981). 

Considérese un modelo W = g(X) que 
involucra una variable aleatoria simple “X”, 
cuyos tres primeros momentos o la función 
de densidad de probabilidad o la función de 
distribución de probabilidad son conocidas. 

En relación con la figura 1, el método de 
estimación probabilística de Rosenblueth 
(EP) se aproxima a la función de distribución 
de probabilidad de la variable “X”, al 
suponer que la masa de la distribución de 
probabilidad de “X” se concentra en dos 
puntos x_ y x,. Usando los dos puntos de 
aproximación, las ubicaciones de x_ y x,, y las 
masas correspondientes p_ y p,, se determinan 
para preservar los tres primeros momentos de 
la variable X. Sin cambiar la naturaleza del 
problema original, es más fácil tratar con la 


: , , (A E 
variable estandarizada X= nn la cual tiene 
media cero y variancia unitaria. Por lo tanto, 
en términos dex”, ,p_yP; 
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f(x) 


Py 


Figura 1. Diagrama esquemático con el método de 


Rosenblueth. 
p+p=1 (20) 
p A +p_ == Q1) 
paó=pt=0"y=1 (22) 
px — px = Y (23) 
En estas ecuaciones, x'= |. Ñ e E 
o] o] 


y yes el coeficiente de asimetría de la variable 
aleatoria “X”. Resolviendo las ecuaciones (20) a 
(23) en forma simultánea se obtiene: 


2 
au=Yy 6 (24) 
2 la 
x=Xx Y (25) 
me 
y (26) 
dd ex 
p_=1-p, (27) 


Cuando la distribución de la variable 


AA 


aleatoria “x” es simétrica, es decir, cuando y = 
0, entonces las ecuaciones (24) a (27) se reducen 
"=1yp_=p,=0.5. Esto implica que para 


ax =xY' 
una variable simétrica aleatoria, los dos puntos 


están localizados a cada lado de la desviación 
estándar de la media, con igual probabilidad 
de masa en los dos puntos. 


A partir de x” y x”, los dos puntos en el 
parámetro original de espacio x_ y x, pueden 
determinarse respectivamente como: 


x_=u-x0 (28) 


x,=4+x/0 (29) 
Basado en x_ y x,, los valores del modelo 
W = g(X) en los dos puntos pueden calcularse 
como: w_=g(x_) y w,=g(x ). 
Entonces, los momentos alrededor del 
origen de W = g(X) de cualquier orden pueden 
estimarse como: 


Elw"]=y0,=p,w'+pw" (30) 

Además del método de la estimación 
de primer orden de la variancia (FOVE), 
el método de la estimación probabilística 
(EP) de Rosenblueth provee una capacidad 
adicional que permite hacer el análisis que 
toma en cuenta la asimetría asociada con la 
función de distribución de probabilidad (DPA) 
de una variable aleatoria. Karmeshu y Lara- 
Rosano (1987) demuestran que el método de 
FOVE es una aproximación de primer orden 
del método de Rosenblueth de estimación 
probabilística. 

Es un caso general donde el modelo 
involucra N variables. El n%"" momento de 
salida W= g(X,, X,,..., X,) alrededor del origen 
puede aproximarse como: 


E(W")-H p(3), 5,,..., By)[w(8, 5... Sy)]” (SD) 


En la que el signo del indicador con 
subíndice $, puede ser sólo + 0 —, representando 
la variable aleatoria, x, teniendo el valor de: 


X¡¿ MF 0 


X=; X_0; 
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Respectivamente, la probabilidad de masa 
en cada uno de los 2* puntos, p(6,, 0,,..., 0,) 


puede aproximarse como: 


N N1A/ N 
pl8, 8, «., 5=[[ 73, 22 50) 62) 
] 


E all y=i+1 


(33) 


Donde p, es el coeficiente de correlación 
entre las variables aleatorias x, y x, El número 
de términos entre la sumatoria de la ecuación 
(32) es 2N, que corresponde al número total de 
posibles combinaciones de + y — para todas las 
N variables aleatorias. 


d) Método de aproximación probabilística de 
Harr. Método del punto de estimación 


Para evitar la intensiva naturaleza compu- 
tacional del método de estimación probabi- 
lística de Rosenblueth cuando el número 
relativamente 
grande, Harr (1989) propuso un método 
probabilístico alterno que reduce el número de 


de variables aleatorias es 


evaluaciones requeridas de 2" a 2N y amplía 
significativamente el método de estimación 
probabilístico a un análisis probabilístico de 
problemas prácticos. 

El método se aplica al segundo momento 
estadístico y es capaz de tomar en cuenta los 
dos primeros momentos (media y variancia) 
de las variables aleatorias involucradas y sus 
correlaciones. Los coeficientes de sesgo de 
las variables correlacionadas se ignoran en el 
método. Las bases teóricas del método de Harr 
se basan en transformaciones ortogonales de 
la matriz de correlación. 

La transformación ortogonal es una herra- 
mienta importante para tratar problemas con 
variables aleatorias correlacionadas. El objetivo 
principal de la transformación es mapear 
las variables correlacionadas de su espacio 


original a un nuevo dominio, en el cual no 
están correlacionadas, por lo que el análisis se 
simplifica significativamente. 

Considérese N variables aleatorias multi- 
variadas X=(X,, X,, 
un vector de valores medios: 


.., X,)' que tienen asociadas 


Ux= (ly, Mo. My)” (34) 
Y una matriz de correlación R(X): 
1 Pp Pa Pin 
] 1 
cam=ro=|P2 7 PB 0 Pal (35) 
Pri Prnz Pms coo 1 


Nótese que la matriz de correlación es 
simétrica; esto es: p, =p, para 14]. 

La transformación ortogonal puede hacerse 
usando la descomposición de los eigen valores 
o también llamada descomposición espectral, 
en la que R(X) se descompone como: 

ROO =C(X”)=VAV! (36) 

Donde V en una matriz de tamaño NxN, que 
contiene N eigen vectores como V = (o, Va.) 
D,), considerando que V, es la 1%" columna de 
los eigen vectores y que A=diag(%, A,..., Ay) 
es una matriz diagonal compuesta de eigen 
valores. 

En términos de los eigen vectores y los 
eigen valores, el vector aleatorio en el espacio 
paramétrico puede 


original expresarse 


como: 
X=u+D'""VA PY (37) 
En la que Y es un vector de N variables 
estandarizadas normales, que tiene a D como 
vector medio y la matriz identidad I, como la 
matriz de covariancia, y D es una matriz diago- 
nal de variancias de N variables aleatorias. 
Las son 
funciones lineales de las variables originales 


variables transformadas Y 


aleatorias, por lo que las variables aleatorias 


Tecnología y 


Ciencias del Agua, octubre-diciembre de 2010 


Ciencias del Agua, octubre-diciembre de 2010 


«o180j0u3aJ 


Marengo ef al., Evaluación de riesgos en proyectos hidráulicos de ingeniería. Incertidumbres y confiabilidad 


X están normalmente distribuidas; entonces 
las variables aleatorias se transforman 
en Y, variables aleatorias estandarizadas 
normales. 

Para un modelo multivariado W = g(X, 
Xp, e en el que se involucran N variables 
aleatorias, el método de Harr selecciona 
los puntos de evaluación localizados en las 
intersecciones de ejes de los N eigen vectores 
localizados en la superficie de una hiperesfera 
multidimensional, que tiene radio /N en el 
eigen espacio y se calcula como: 

X =4+4/ND"V, ¡=1 2... N (88) 

En la que X, representa al vector de 
coordenadas de las N variables aleatorias en el 
espacio paramétrico, que corresponde al ¡%”" 
eigen vector V, u = (4, U,..., 1, )', un vector de 
medias de N variables aleatorias X. 

Basado en los 2N puntos determinados por 
la ecuación (38), se pueden calcular los valores 
de la función. Entonces el momento n%'” del 
modelo de salida W alrededor del origen puede 
calcularse con las siguientes ecuaciones: 


m m 
a wW_ +w, Mx, "Mx. 
wi = a Le uds (x;-) , para 


(39) 
SL Da Ni mel, Des 


N, —m 
A. 
E[w” =u4"(W)= Za» i, para m=1 2,... (40) 


Alternativamente, la transformación orto- 
gonal puede hacerse con la matriz de 
covariancias. 


e) Método del valor medio del primer orden 
del segundo momento (MFOSM) 


En el método del primer orden del segundo 
momento estadístico, la función de comporta- 
miento, W(X), definida con base en las funcio- 
nes de carga y resistencia, se expanden en series 


de Taylor en un punto dado de referencia. A 
partir del segundo término de la serie se trunca, 
resultando en una aproximación que requiere 
del conocimiento de los dos primeros momen- 
tos estadísticos de las variables aleatorias. Esta 
considerablemente 
la practicidad del método de primer orden, 


simplificación aumenta 
porque en muchos problemas reales es muy 
difícil, si no imposible, encontrar la DPA de las 
variables que intervienen en los problemas; sin 
embargo, es muy fácil obtener los dos primeros 
momentos estadísticos de dichas variables. 
Los procedimientos para estimar el método de 
primer orden del segundo momento estadístico 
se muestran con detalle en la literatura 
especializada que describe el método FOVE 
para el análisis de incertidumbres. 

Una vez que se han estimado la media y la 
desviación estándar de W(X), la confiabilidad 
se expresa como: 


Mo 


Burosm 3 TcEDs 7 COS 


(41) 


Dónde u y C(X) son los vectores de medias 
y matriz de covariancia de las variables 
aleatorias X, s=V.W(u) 
es el vector columna de los coeficientes de 


respectivamente; 


sensibilidad de cada elemento que es evaluado 
enX =p. 


f) Aproximación de primer orden de la función 
de comportamiento en el punto de diseño, 
método avanzado del primer orden del segun- 
do momento (AFOSM), según Hasofer-Lind 


La principal ventaja del método AFOSM es 
tratar de mitigar las deficiencias asociadas 
con el método del primero orden del segundo 
momento estadístico (MFOSM) anterior, pero 
conservando la simplicidad de la aproximación 
del primer orden; es decir, calculando sólo 
los dos primeros momentos estadísticos de 
las variables estudiadas. La diferencia entre 
el MFOSM y el AFOSM es que el punto de 
expansión de primer orden en este último 
está localizado en la superficie de falla y el 
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primero en el punto medio de las variables 
estandarizadas, por lo que los investigadores 
que estudian estos temas recientemente 
prefieren el método avanzado en lugar del 
método de primer orden. 

Considerando el punto de expansión en X, 
=X: 

WOO=Y 55 (xj) =sH(X-Xp) (42) 

En la que s, = (s,., S,...Sy») es un vector de 
sensibilidad de los coeficientes de la función de 
comportamiento W(X) evaluado en el punto de 
expansión X,. 

El valor de S es: 


a oW(X) 
! 0x; Sica 


Entonces el valor medio y la variancia de la 


(43) 


función W(X) pueden expresarse como: 


y =s:(1—X:) (44) 


oy =s COOs. (45) 
En la que C(X) es la matriz de covariancias 
de las variables aleatorias; si las variables no 
están correlacionadas: 
N 
Ow=)Y 550, (46) 
La desviación estándar de la función de 


comportamiento W(X) puede expresarse en 
términos de las derivadas direccionales como: 


N 


Donde a... es la derivada direccional para la 
1 
variable 1%" en el punto x.: 


cy N (48) 


En forma matricial: 


DYVW(X,) 


OL ¡+ ETE NAREO Y 49 
ID'PVW(X) a 


El índice de confiabilidad se estima enton- 
ces como: 


N 

—Hyw_ 2 Si (1;—Xz) 
= N 

Sw Dd 0rS»0; 


(50) 


Una vez que B se ha calculado, la 


AFOSM 
confiabilidad puede estimarse como p_ = 
D(Barosu); la sensibilidad de las variables se 
calcula como: 


E ON 
pe! 


Obs (51) 


Ve Barosm = Vr |x 


Ecuación que muestra que ol, es la relación 


del cambio en f debido a un cambio en la 


AFOSM 
desviación estándar de la variable X, en X = 
X.. Entonces la relación entre V .B y VB puede 
expresarse como: 

Ve Barosm=D'"?Vy Barosm=-D "Wa. (52) 

También puede demostrarse que el 
coeficiente de sensibilidad de la confiabilidad 
o probabilidad de falla con respecto a cada 
variable estocástica puede ser expresada 
como: 


e =D, rl e 6 
X: 


9p.1_ Urb Barom. ¡12 (54) 
OX, O; E dl A Ed 
En forma matricial: 
Vo ps= 0 (Barosm)V: Barosm 
(55) 


= 0 (Barosm)D'Pa.- 


El coeficiente de sensibilidad muestra 
la importancia relativa de cada variable 
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estocástica en la confiabilidad o probabilidad 
de falla. 

El algoritmo del método avanzado del segun- 
do momento para variables independientes 
según Hasofer y Lind (AFOSM) es el 
siguiente: 

Si X son variables normales independientes, 
la estandarización se reduce a Z variables con 
media cero y matriz de covariancia Í, siendo I 
una matriz de identidad NxN. Hasofer propuso 
la ecuación recursiva: 


ze =-|-a, zip Jo WE 10) 

(1+1) (1)4(1) 11) VW(Z¿,) (1) (56) 
PEL. Z ds 
Donde: 


r, esla iteración ¡$ 
a es el vector unitario de la superficie de 
falla. 


En el espacio, x se escribe: 


o -1] sc) WX) 


X= 1+Dsp 
(+0 (o) 
S(» DSs(» 


(y (57) 
r=L. Bis 


Basado en esta ecuación, el algoritmo de 
Hasofer-Lind con el método AFOSM queda 
como sigue: 


1. Elegir una solución de prueba X,,.. 
Calcular W(X,)) y el correspondiente vector 


Se 


3. Revisar el punto de solución X,,.,, 


de 
acuerdo con la ecuación (57). 


4. Verificar si Xo) y Xy son lo suficientemente 


cercanos; en caso afirmativo, calcular B 20 


con la ecuación: 


N 
Mw _ A = X;») 


= (58) 
Sw DD .0ys):0; 


Barosm 3 


Y el índice de confiabilidad P.= p(B 


AF A) 


Si X,, * X 


(Hay entonces repetir 2; en caso 


contrario: 


5. Calcular el índice de confiabilidad con 
respecto a los cambios de las variables 
estocásticas. 


Debido a la naturaleza no lineal de la 
optimización, el algoritmo anterior no converge 
necesariamente en el verdadero punto de 
diseño asociado con el mínimo índice de 
confiabilidad, por lo que Tang propuso un 
método que lo vuelve convergente a dicho 
punto de diseño. 


g) Método avanzado del primer orden del 
segundo momento estadístico según Tang 


método 
del primer orden del segundo momento 


La formulación matemática del 


estadístico se muestra con detalle en el libro 
de Tang (1980); sin embargo, a continuación 
se presentan las bases de la formulación 
matemática y un resumen de su aplicación. 

Los conceptos relativos a confiabilidad 
pueden limitarse a una formulación basada en 
el primero y segundo momentos estadísticos 
de las variables aleatorias que intervienen en 
el problema de la formulación del segundo 
de ellos (Cornell, 1969; Ang y Cornell, 1974). 

Con el enfoque del segundo momento, la 
confiabilidad puede medirse completamente 
en función del primer y segundo momento de 
las variables de diseño. 

Tomando en cuenta que la definición de 
margen de seguridad es M = X - Y, el “estado 
de seguridad” se define para M>0 y el “estado 
de falla” para M<0. La frontera que separa los 
estados de falla y seguridad queda establecida 
para M = 0. Si se consideran las variables 
reducidas: 


X'= y Va 


x Sy 


(59) 


El espacio de estas variables reducidas se 
muestra en la figura 2. También en términos de 
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Estado de falla 
M<0 


d 


Estado seguro 
M>0 


Figura 2. Espacio de las variables reducidas X' y Y.. 


las variables reducidas, la ecuación límite M = 

0 viene a ser: 
0,X=0,Y'+4,1,=0 (60) 

que es la línea recta mostrada en la figura 2. La 


distancia desde la línea de falla al origen 0 está 
dada por: 


(61) 


Esta expresión se cumple solamente si las 
funciones de distribución de probabilidad de 
las variables analizadas son normales y no 
están correlacionadas entre sí. 

Lo anterior sucede porque sólo con la 
distribución normal, en el caso de variables 
no correlacionadas, se cumple que la media 
de la suma de las variables que intervienen 
en el problema es la suma de las medias de 
la función de comportamiento; lo mismo se 
aplica al hecho de que en este caso la variancia 
de la función de comportamiento es la suma 
de las variancias de las variables estudiadas. 

Debido a lo anterior, en el caso de que 
se tengan variables con distribuciones de 
probabilidad diferentes a la normal, se debe 
evaluar la probabilidad que corresponde con 
base en la distribución normal equivalente, 
cuya aplicación se muestra en diversos libros 
de estadística. 


Ya 
interpretarse como una 
confiabilidad del 
evaluación de la ecuación (61) en el origen de 
las variables reducidas (M' = 0) al tratarse con 


que la distancia mínima puede 
medida de la 


sistema, entonces la 


distribuciones normales puede hacerse de la 
siguiente manera si: 


M=X-Y (62) 

> My = Hx —Hy (63) 
0%,=02-0? (64) 
M'=(M-um)/0u (65) 


En el origen de M(M' = 0), para M=0 
(distancia entre el origen y el estado de falla 
según la figura 2) conduce a que: 


Pr = Ey (MO = Em (0)= Fun (Hu! 90m) (66) 
= 0 (4 1/0) =1-6(Um 15) 


El índice de confiabilidad $ es esta distan- 
cia mínima al origen (u,,/0,,), o sea que: 


B=H1/0u 


P¿=1-P¿=1-1+4 (P)=0 (PB) (67) 

Esta demostración hecha para dos variables 
(X y Y) en el plano puede generalizarse a tres 
dimensiones O al espacio de n variables no 
lineales, al encontrarse un plano tangente a 
la superficie de falla y la distancia de éste al 
origen. 

Los resultados relevantes de la formulación 
del método avanzado de primer orden de los 
segundos momentos pueden resumirse como 
sigue: 


e El punto más probable de falla se puede 
calcular con la ecuación: 


X/=-0p (68) 
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en la cual a, son las direcciones de los cosenos 


08 
oX; 


directores: 


(69) 


donde las derivadas son evaluadas en (X;*, 
DE y 


X= O,¡X¡ + Mi = My; 0%; Ori B (70) 
e La solución de la ecuación límite de estado 
permite obtener f. 


Los resultados resumidos anteriormente 
permiten plantear el siguiente algoritmo para 
aplicar el método avanzado de primer orden 
de los segundos momentos estadísticos, y que 
es el empleado en el cuerpo del trabajo: 


1. Definir una función de comportamiento 
“o(X)” con las variables estadísticas que se 
consideren adecuadas y significativas en el 
problema analizado. 

2. Suponer un punto inicial de falla X'*; ¡ =1, 
2, ...,..., M y Obtener: 


X'= Et — Pas (71) 


Xi 


3. Estimar (0g/0X'*) y a en los puntos x/. 

4. Calcular x/ con: X =H,,—0o,.B. 

5. Sustituir los valores estimados de X/ en la 
función de comportamiento y encontrar Pf. 

6. Conel valor de B obtenido reevaluar X'* =- 
op. 

7. Repetir los pasos tres a seis hasta que 
se obtenga la convergencia deseada que 
sucede cuando B,=B.,,. 

Debe señalarse que con la aplicación de este 
método no es posible estimar el valor medio de 
la función en conjunto debido a la naturaleza 
no lineal de la aplicación del método. 


Aplicación 


Para mostrar la aplicación de los métodos 
indicados se busca determinar la confiabilidad 
para poder compararlas entre sí. 
utilizando la 
ecuación de Manning, con las siguientes 
condiciones (Mays y Tung, 1992): 


Una alcantarilla trabaja, 


Qp=0.4631*D*%505 (72) 


Donde: 


Q = gasto en ft /s. 

n = coeficiente de rugosidad en ft'/*. 
D = diámetro en ft. 

S = pendiente del tubo. 


Las propiedades estadísticas de las varia- 
bles se muestran en el cuadro 1. 

La confiabilidad debe calcularse para cono- 
cer las condiciones en las que la alcantarilla 
descargue 35 ft” /s. 


Método de la integración directa 
Para la alcantarilla analizada, con variables 
aleatorias n, D y S, se tienen las siguientes 


propiedades estadísticas (cuadro 1). 


La función de resistencia es: 


W(n, D, S)=0.0463 1 D?50 (73) 
Y la función de carga es L = 35 
>W(n, D, S)=In(R)-In(L) (74) 


Cuadro 1. Propiedades estadísticas para el Método de 


Integración Directa. 
Parámetro Media Coeficiente de variación 
n 0.015 0.05 
D 3.00 0.02 
S 0.005 0.05 
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Considerando que las variables n, D y S son 
independientes y lognormales, la función de 
comportamiento puede escribirse: 


W!(n, D, S) 

=[In(0.463) - In(1) +2.67 In(D) +0.5In(S) |-In(35) 
W(n, D, S) 

= —In(1) +2.671n(D)+0.5 In(S)-— 4.43319 


p.=p[W(n, D, S)>0] (75) 


n D y S son variables lognormales 
independientes; debido a que In(n), In(D) y 
In(S) son variables aleatorias independientes, 
entonces la función de comportamiento W(n, 
D y S) es una función lineal de variables 
aleatorias independientes; entonces W(n, D y 
S) es una función lineal de variables aleatorias 
normales independientes, cuya media, de 
acuerdo con la propiedad reproductiva de 
variables aleatorias es: 


(76) 
=0.15853 
Y variancia: 
Q%=-Q% ,+(2.67)Varln(D)+(0.5)*Varln(S) n 
77 


=0.005977 


La confiabilidad se puede obtener como: 


_ al Mw )]_,(_0-15853 
p=p00>0=0[ 3%) o ] 


+/0.0059765 
=0(2.050607) =0.97818 
Pr =1-0.97818 = 0.02182 
Método de Monte Carlo 
Una vez que se establece la función de 


comportamiento respectiva (ecuación (73)), se 
genera una muestra significativa de números 


aleatorios de cada variable y se estima cuántas 
veces excede el valor cero; el riesgo de falla 
se evalúa como la relación entre el número 
de valores positivos o negativos de F, (0) con 
respecto al total de elementos generados. 

Se adoptaron las características estadísticas 
de media y desviación estándar para cada 
variable, con distribuciones de probabilidad 
normal para cada una de ellas, y se generaron 
3 000, 4 000, 5 000 y 6 000 números aleatorios, 
con una distribución de probabilidad uniforme 
en un rango U(0,1), obteniendo como resultado 
para cuatro corridas diferentes los resultados 
mostrados en el cuadro 2. 

Para el análisis se adoptó como resultado 
el obtenido para 6 000 números aleatorios, que 
dio como resultado una probabilidad de falla: 


p,=0.02011 


Como se señaló anteriormente, el método 
de Monte Carlo es sensible al tamaño de 
la muestra y aún con un gran número de 
simulaciones, 


en cada caso se obtienen 


resultados diferentes. 
Método de Mellin 


En el ejemplo de aplicación considere que 
los tres parámetros del modelo son variables 
aleatorias independientes con las propiedades 
indicadas en el cuadro 3. 

De acuerdo con la transformación de Mellin, 
el gasto M¿(s) puede obtenerse como: 


M/(s) = 0.463 M,(=s+2)Mp(0.55+5) (78) 


Para el coeficiente de rugosidad n, teniendo 
distribución uniforme, se obtiene: 


3 +2 


Ey A 


=s+2_ 


ny n 


M,(=s+2) = 


Para el diámetro de la alcantarilla con una 
distribución triangular se obtiene: 
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Cuadro 2. Resultados de probabilidad de falla obtenidos con el método de MonteCarlo con muestras 


aleatorias de 3 000, 4 000, 5 000 y 6 000 números generados con distribución de probabilidad uniforme. 


ena Corrida 1 Corrida 2 Corrida 3 Corrida 4 Promedio 
3 000 58 /3 000=0.01933 62/3 000=0.02066 68/3 000=0.02266 57 /3 000=0.01933 0.02050 
4 000 88 / 4 000=0.02200 83 / 4 000=0.02075 80/4 000=0.02000 85 /4 000=0.02125 0.02100 
5 000 107/5 000=0.02140 109/5 000=0.02180 101/5 000=0.02020 106/5 000=0.02120 0.021150 
6000 125 /6 000=0.02083 126/6 000=0.02010 122 /6 000=0.02033 115/6 000=0.01917 0.02011 
Cuadro 3. Propiedades estadísticas para estimar la confiabilidad con la transformación de Mellin. 
Parámetro Distribución 
n Distribución uniforme con límite inferior de 0.0137 y superior de 0.0163 
D Distribución triangular con límite inferior de 2.853, media 3.0 y superior de 3.147 
Distribución uniforme con límite inferior de 0.00457 y superior de 0.00543 


a 
(d, -d,12.67s-1.671(2.675-0.67) 
di ld ¿895-1672 de e5Lon d, (4 2% -1.67_ ae =1.67) 


d,—d d, —d 


m qa 


Mp(2.675-1.67)= 


m 


(80) 


Para la pendiente de la alcantarilla S con 
distribución uniforme se obtiene: 


(055+05) _ (, y (05s+0,5) 
Ms(0.5s +0.5)= 0) (s,) 
(0.55 +0.5)(s, =5,) 


Basado en el anexo Ll, la transformada 
de Mellin para cada parámetro del modelo 
estocástico se puede expresar como se resume 
en el cuadro 4. 

El valor medio de la capacidad de la 
alcantarilla puede determinarse como: 


E(Q)=Mp(s=2) 
=0.0463M,(0)Mp(3.67)Ms (1.50) 
= 0.463 (66.834) (18.806)(0.0707) =41.137 £té/s 


El segundo momento sobre la capacidad 
original de la alcantarilla es: 


E(Q?)=Mpg(s =3) 

= (0.0463)”[M, (DIM p(6.34)M5(2.00) 
= (0.463)*(4 478.08)(354.681)(0.005) 
=1702.40 fté/s 


La variancia de la capacidad de la alcanta- 
rilla puede determinarse como: 


Var(Q) = E(Q?) - E*(Q) =1 702.40 - (41.137)? 
=10.1472 (£é/s)? 


Siendo la desviación estándar: 


O9= ÍVar(Q) = 410.147 =3.185 ft%/s 


La media y desviación estándar de la función 
de comportamiento W son, respectivamente: 


yy =H9-35=41.14-35=6.14 ft/s 
Oy =07=3.185 ft/s 


6.14 


=— = 1.928 
3.185 
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Cuadro 4. Parámetros de cálculo con la transformada 


de Mellin. 
s=2 5=3 
0.463! 0.463 0.2144 
M,(+s+2) M,(0) = 66.834 M, (+) = 4 478.080 


M,.67-1.67) | M,(3.67) = 18.806 | M,(6.34) = 354.681 


M,(0.5s + 0.5) | M.(1.50) =0.0707 | M6¿(2.00) = 0.005 


p¿=0 (1.92476) =0.973197 
La probabilidad de falla es: 
p;=1-ps=0.0268 


Técnica aproximación: método de 
estimación de primer orden de la 
variancia (FOVE) 


no están 
correlacionas entre sí, la expansión en series 


Suponiendo que las variables 


de Taylor de la función de comportamiento 
alrededor de 1, u, = 0.015, D, 1, = 3.00 y S, hu, 
= 0.005, es: 


W=0.463 (0.015)7—(3)2%(0.005)%5+4 ES (n—0.015) 
n 


9Q -PaQ 
+ lo 0 Els 0.0005) 


0D 


W= 41.01+| 0.463 (-D(0.015) (3.0 y%(0.005)% ] 
(n-0.015) 
+| 0.463(2.67)(0.015)(3.0)"*(0.005)% |(D-3.0) 


+[0.463(0.5)(0.015) (3.0)?%(0.005 ) 05] (S-0.005) 


W= 41.01 - 2 733.99 (n — 0.015) + 36.50 (D-3.0) 
+ 4 100.99 (5 —0.005 ) 


Tomando en cuenta la ecuación (17), la 
aproximación de la media de la capacidad del 
flujo es: 

y = 40.96 £t/s 


De acuerdo con la ecuación (18), la 
aproximación de la variancia de la capacidad 
del flujo es: 


Oy = (2733.99)"Var (n) + (36.50) Var (D) 
+ (4 100.99) Var(S) 

-2(2 733.99)(36.50)Cov(n, D) 

-2(2 733.99)(4 100.99)Covín, S) 
+2(36.50)(4 100.99)Cov(D, S) 


Esta expresión se reduce a: 


Oy = (2733.99)? Var(n) + (36.50)? Var(D) 

+ (4 100.99)? Var(S) 

Debido a que Covín, S) =Cov (D, S) =Covín, 
D) =0. 


Como las desviaciones estándar de la 
rugosidad, diámetro del tubo y pendiente son: 


o, = (0.05)(0.015) = 0.00075 
Op =(0.02)(3.0) = 0.06 
Os = (0.05)(0.005) = 0.00025 


La variancia de la capacidad de la alcanta- 
rilla se calcula como: 


oy =(2733.99)*(5.625x107)+(36.50)*(3.6x107) 
+ (4 100.99)*(6.25x10*) 

-2(2 733.99)(36.50(-3.375x107) 

07 =2.057+2.19+1.03?+6.74 =10.05172 (0 


La desviación estándar de la capacidad del 
flujo es: 


Oyy= 110.05172 =3.17045 ft/s 
La media del gasto es: 
yy =40.96-35=5.96 ft%s 


Sy = v10.05172=3.17045 té/s 


5.96 
=0 | 2% | =0(1.8798)=0.9699 
ds (517055) ¡id 


La probabilidad de falla es: 


pp=1-0.9699 = 0.0301 
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Método de Rosenblueth 


En cuanto al problema de la alcantarilla 
analizada, supóngase que los tres parámetros 
de la ecuación de Manning son variables 
aleatorias. 

Considerando la fórmula de Manning en la 
que: 


W=0Q=0.463nD?%503_35 


Con tres variables aleatorias hay un total de 
2? =8 posibles puntos a ser considerados con el 
método de Rosenblueth. Debido a que hay tres 
variables simétricas, el coeficiente de asimetría 
es igual a cero, por lo que de acuerdo con las 
ecuaciones (24) a (27): 


, , 


n=1, 


D'=D,=S'=S/=1 


Y los valores correspondientes del coefi- 
ciente de rugosidad, diámetro de la alcantarilla 
y pendiente son: 


n, =p, +10, =0.015+(1)(0.00075)=0.01575 
n_=u4,-10,=0.015- (1)(0.00075)=0.01425 
D,=1p+ Do, =3.0+(1)(0.06)=3.06 ft 
D =up-D'op=3.0-(1)(0.06)=2.94 ft 
S, =15+5S/05=0.005+ (1)(0.00025)=0.00525 
S_=45-S'0,=0.005- (1)(0.00025)=0.00475 


Sustituyendo los valores de 1”, 1”, D', D”, S” 
y S” en la ecuación de Manning, se obtiene, por 
ejemplo: 


Qi= 0.463(n,J UDS, 95 


=0.463(0.01575)(3.060.00525)= 42.19 f£é/s 


Similarmente, los valores de la capacidad 
de la alcantarilla para los otros siete puntos 
están dados en el cuadro 5. 

Debido a que el coeficiente de rugosidad y 
el diámetro de la alcantarilla son simétricos, 
las variables correlacionadas de las variables 
aleatorias las 


pueden determinarse con 


ecuaciones (32) y (33). 


Puig "Po =(1+P,p+ Pas +Pos)/8 
= (1+0+0+0)/8=0.125 
Pi Pos =(l+Pp7 Pas”Pps)/8 
=(1-0-0-0)/8=0.125 
Pues Ps =(1-Pp+ Pas—Pns)/8 
= (14+0+0-0)/8=0.125 
Prr =P4- =(1-Prp7 Pas +Pps)/8 
=(1+0+0+0)/8=0.125 


Los valores de probabilidad en masa 
tabulados se muestran en la última columna. El 
momento de orden m%'"" alrededor del origen 
para determinar la capacidad del flujo puede 
determinarse con la ecuación (31), como se 
muestra en el cuadro 6. 

Del cuadro 6, 1, = E(Q) = 41.22 fé /s y E(Q?) 
=1 715.994 (ft /s)?. 

Entonces la variancia de la capacidad del 
flujo de la alcantarilla puede estimarse como: 


Cuadro 5. Capacidad de la alcantarilla con el método de Rosenblueth. 


Punto n D Ss Q (0) 
1 + 42.19 0.03125 
2 = 40.14 0.03125 
3 - 37.92 0.21875 
4 = = 36.07 0.21875 
5 - 46.64 0.21875 
6 = = 44.36 0.21875 
7 = - 41.19 0.3125 
8 = = = 39.81 0.3125 
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var(W)= E(Q?) - (10) =1 702.419- (41.219)? 
=10.022 (fe Js) 


Por lo tanto, la desviación estándar de 
la capacidad de flujo de la alcantarilla es 
y10.022=3.166 ft/s. 

Comparando con los resultados anteriores, 
se Observa que el método de estimación proba- 
bilística de Rosenblueth arroja valores más 
altos de media y variancia que los obtenidos 
con el método de estimación de primer orden 
de la variancia (FOVE). 

Para estimar la confiabilidad: 


lyy=41.139 (£té/s) 
0?,=10.022 (£6/s) 
> ly =41.139-35=6.22 (fs); 
Oy =3.166 (16%) 


E O E 


wW 


al A _ 
lares 0 (1.939) = 0.97381 


La probabilidad de falla es: 


Pp=1-0.97381=0.02619 


Método de aproximación probabilística 
de Harr. Método del punto de estimación 


De acuerdo con la ecuación (38), las coorde- 
nadas de los 2x 3=1 puntos de intersección 
que corresponden a los tres eigen vectores y la 
hiperesfera con radio /3 pueden determinarse 
como: 


Us oc, 0 0 0.015 
Xw=|up|+43|0 op 0|V=| 3.0 
Us 0. 0. 06s 0.005 
0.00075 0 0 
+43| 0 0.06 o  |v, 
0 O  0.00025 
parai=1,2,3 
0.7071 
Dónde V, es un vector con valor Y,=| 0.7071 
0.7071 


ya que las variables no están correlacionadas. 

Las coordenadas resultantes de los seis 
puntos de intersección de estas ecuaciones 
se enlistan en la columna 2 del cuadro 7. 
Sustituyendo valores de x (columna 2) 
en la fórmula de Manning se calculan los 
correspondientes al gasto (columna 3); los 
valores de Q* se presentan en la columna 4, 
calculando el segundo momento respecto al 
origen. 


Cuadro 6. Cálculos de los primeros dos momentos alrededor del origen. 


Punto Q P Q*P Q*P Q>P 
(1) (2) (3) (4) (5) (6) 
1 42.19 0.125 5.274 1779.9961 222.5 
2 40.14 0.125 5.018 1611.2196 201.42 
3 37.92 0.125 4.740 1 437.9264 179.741 
4 36.07 0.125 4.509 1 301.0449 162.631 
5 46.64 0.125 5.830 2 175.2896 271.911 
6 44.36 0.125 5.545 1967.8096 245.976 
7 41.91 0.125 5.239 1756.4481 219.556 
8 39.87 0.125 4.984 1 589.6169 198.702 

Suma 1 41.139 1702.419 
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Se obtienen las columnas 3 y 4, y los valores pp=1- ps =0.02438 
promedio de Q y Q? a lo largo de cada eigen 
vector son calculados y listados en las colum- Método del valor medio del primer 
nas 5 y 6, respectivamente. orden del segundo momento (MFOSM) 
De acuerdo con la ecuación (40), con m=1 se 
puede calcular el valor medio de la capacidad Los procedimientos para estimar el método 
de la alcantarilla como: de primer orden del segundo momento esta- 
dístico se muestran con detalle en la literatura 
_MQ;+2,0,+1Q3 41.143 10/s especializada que describe el método FOVE 
ld MAA +A para el análisis de incertidumbres. 
Una vez que se han estimado la media y la 
El segundo momento respecto al origen se desviación estándar de W(X), la confiabilidad 
calcula como: se expresa como: 


=> —2 2 
E(Q?) e MQí +10) +20; =1 702.432 (1e/sy Hi 


A+ +A Puros ¡Isis 
La variancia del flujo en la alcantarilla es: Con el método FOVE se estimó que: 


ES 2 2_ 
var(W)= E(Q5)-(1,)=1702.432 e La media del gasto y la desviación estándar 
- (41.143?) =9.6612 son: 


llyy= 41.143 -35= 6.143 ft/s; Hi =40.96-35=5.96 (£t/5) 
var (W)=9.6612 (Ets); 


Oy =+10.05172 =3.17045 (£t/s) 
Oy = 49.6612 = 3.10855 £'/s 


5.96 
.14 = 0 | ——— |=0 (1.8798) =0.9699 
Ber == 30055 1.9763 Ps (537085) , ) 
p.=0(1.9763) =0.97562 e La probabilidad de falla es: 
La probabilidad de falla es: p¿=1-0.9699 = 0.0301 


Cuadro 7. Evaluación de los momentos estadísticos con el método de Harr. 


Punto X= (n, D, S) Q Q? Q_ (0 
(1) (2) (3) (4) (6) 
E > | (oowosj3oz35 0005 | «6606 | 27218 [| | 
E 3 | Cos53ooo0s8m) | «26 | 17606 | | 
3 3- (0.015, 3.00, 0.00469) 39.718 WTI 
8 
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Método avanzado del primer orden del 
segundo momento (AFOSM) según 
Hasofer-Lind 


Tomando en cuenta el algoritmo de Hasofer- 
Lind con el método AFOSM, la ecuación de 
comportamiento se puede expresar como: 


Wín, D, S)=Q7-Q¿=35-0.463nD?*%505 
El algoritmo que se sigue es: 


1. Elegir una solución de prueba X,,, con los 


valores medios iniciales: 
X=(l: Mp Hs) =(0.015,.-3.00, 0.005)" 


2. Calcular W(X;) y el correspondiente vector 
Ss 
(1) 


Win, D, S)=-6.010 


E area! ; S,=2734; Sp =0.00365; 
X=x* 


Sg =0.00004101 


3. Revisar el punto de solución de acuerdo 
con la ecuación de comportamiento: 


Wí(n, D, S)=-6.010 
y 
dif=|X,-X,||=[(0.01592 -0.015)?+(2.921-3.0)? 
+(0.004847 -0.005)*]””=0.07857 
4. Verificar si X,,, y X,,,, son lo suficientemente 


cercanos; en caso afirmativo, calcular f 
con la ecuación: 


My _ E 7; — 59) 
Sw E Si “O; 


Y el índice de confiabilidad P, = p(B 


AFOSM 


Bar OSM 


dea 


Si X,, ES Xy 
Pl 


y entonces repetir 2; en caso 


5. Calcular el índice de confiabilidad con 
respecto a los cambios de las variables 
estocásticas. 


El cálculo se muestra en el cuadro 8. 
En este caso, la media y desviación estándar 
del gasto son: 


¿(M¿— 2%) =5.536 fs; O ,, 


My = Ys 
DI 


29 [5536 
Ps 2.691 


s,*0¡=2.691 ft /s 


)- 0 (2.057) =0.9802 


La probabilidad de falla es: 


p=1-0.9802=0.01983 


La sensibilidad del método se muestra en el 
cuadro 9. 

By 
ox; ox; 
sensibilidad del indlies de 
confiabilidad por el cambio de una vez de la 


muestran las 


desviación estándar, mientras que muestra que 


%B y op, corresponden al cambio unitario de 
X; 


ox, > 
las variables aleatorias en el espacio original. 
La sensibilidad de PB y p, asociadas con 
el coeficiente de Manning son negativas; 
mientras que las asociadas con el diámetro 
y la pendiente son positivas; esto indica que 
un incremento en el coeficiente de Manning 
resultaría en un decremento de $ y p,, mientras 
que un incremento en la pendiente y en el 
tamaño del tubo incrementarán $ y p.. 
Hidráulicamente esto se explica, ya que 
un incremento en la rugosidad provoca 
necesariamente un decremento en el gasto, y 
un incremento en el diámetro o la pendiente 
significa una en la 


mayor capacidad 


conducción. 
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Cuadro 8. Índices de confiabilidad del método AFOSM según Hasofer. 


Iteración Variable A Sw Oli) ón 
n 0.015 2734 0.6468 0.01592 
D 3.00 -36.50 0.6907 2.921 
r=1 S 0.005 -4101 0.3234 0.004847 
dif=0.07857 W=-6.010 PB=0.0 
n 0.01592 2226 0.6138 0.01595 
D 2.921 -32.39 0.7144 2.912 
r=2 S 0.004847 -3656 0.3360 0.004827 
dif =0.004847 W=-0.04421 [PB =1.896 
n 0.0155 2195 0.6118 0.01594 
D 2.912 -32.09 0.7157 2.912 
r=3 S 0.004827 -3625 0.3369 0.004827 
dif =0.000001919 W=-0.00251 PB =2.056 
n 0.0154 2195 0.6119 0.01594 
D 2.912 -32.10 0.7157 2.912 
r=4 S 0.004827 -3626 0.3369 0.004827 
dif =0.000003721 W= -0.0000002544 PB =2.057 
Cuadro 9. Sensibilidad del método AFOSM según Hasofer. 
Variable Xi) Ol 9B 9, 9 9Ps x0B xp, 
0 0) 6) Ox; Ox; 9x; 9x; oPx, p,0x, 
(4) (5) (6) (7) (8) (9) 
n 0.01594 0.6119 0.6119 -0.02942 -815.8 -39.22 6.323 0.638 
D 2.912 0.7157 0.7157 0.03441 11.9 0.57 16.890 1.703 
S 0.00483 0.3369 0.3369 0.01619 1347.0 64.78 3.161 0.319 
Método avanzado del primer orden del C11=0.463; C21=-1.0; C3=2.67; C4=0.5 
segundo momento (AFOSM) según Tang 
Las derivadas parciales del gasto son: 
La función de comportamiento del problema 
Ss d 1 1 
sI: 0.463 DM) =C11" A 
n n n 


_ 15267 (05 
Qp=0.463n"D*”S a CUAD IS Ap 02) 


Wí(n, D, S)=35-0.463n"D*%5% 90 
¿y 0-463175"2(2.67D'%) 


Con los parámetros estadísticos: =C11*n 542.67 D%) 


=C11*ns 4c3p (3D 
u,=0.015; 1y=3.00; 1s=0.005 


5, =0.00075; 0 =0.06; 0 5=0.00025 20 —0.4681"D%B(0.557") 


9S 


Si se adoptan como constantes: =P ADA 
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Y la función de comportamiento con el 
margen de seguridad es: 


W(n, D, S)=35-0.463n1D?9505 

La primera iteración se hace considerando 
los valores medios de la variable: 

W(n, D, S) =35-41.00988 = -6.00988 >>0 

Las parciales son: 


1 


9Q 8/3 1/2 
= =(.463(3 0.005 =- 
(37% pel TO 


n 


37) =-2 733.992 


me =0.463 (0.015) (0.005)? (2.67(3)P)=36.499 


e =0.463 (0.015) (3)% (0.5(0.005) 2) 


=4 100.98787 


ro ,0,7 00 2.1/2 
sumM=[( Sn de Has? ] 4 928.909 


»- 0Q/) _ 05547, q CRÍA) 


1” SsuM SsuM 


=0).00741; 0'= Lea = 0.83203; 


Wín, D, S)=35-Q,=35-94% pesgun 
n 


El cálculo se resume en el cuadro 10. 
La ecuación de falla es: 


0.463 


Wín, D, S)=35 
só ) (0.0162 + 0.003512 B) 


(2.99875 — 0.00042f)/* (0.0042 -0.0022B)P=0 


Que se resuelve para $ = 3.39211. 
Para ejemplo en particular, la solución con 
el algoritmo de Tang es: 


1. =0.0162; 01, =-0.4649 
D'=2.9986; 0.,=0.0067 
S'=0.0042; a3= 0.8853 


El coeficiente de confiabilidad es $ =3.39211: 


> P,=0.000348 
P¿=0.999654 


Discusión de resultados 


El resumen de los métodos de estimación de 
riegos se muestra en el cuadro 11. 

Se pueden señalar las siguientes observa- 
ciones: 


1. Se puede adoptar como método exacto 
el obtenido con el de integración directa; 
sin embargo, hay que considerar que 
se supuso que los resultados obtenidos 
obedecen al hecho de adoptar que las 
variables estudiadas tienen una función de 
distribución lognormal; en casos analiza- 
dos con incertidumbres más grandes que 
las estudiadas, este método puede no tener 
aplicación. 


2. Elmétodo de Monte Carlo arroja resultados 


parecidos al de integración directa, con una 
diferencia AP, = 0.07852 (7.852%); puede 
usarse como referencia para la estimación 
del riesgo, aunque en los casos analizados 
que poseen variables con distribuciones 


Cuadro 10. Puntos de falla del método AFOSM según Tang. 


Punt t N t fall 
Variable unto supuesto $ 2Q/ox, E uevo punto de alla 
de falla vi z bi a IO A | 
n 0.015 0.00075 -2 733.992 0.5547 0.015— (-0.5547)(0.00075)B 
D 3.00 0.06 36.499 0.00741 3.00-(0.00741)(0.06)B 
S 0.005 0.00025 4 100.98787 0.83205 0.005-—(0.83205)(0.00025)B 
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de probabilidad diferentes a la normal, el 
análisis presente una mayor complejidad 
en su evaluación. 


3. Si el resultado exacto se adopta como el 


obtenido con el método de integración 
directa, el que más se parece es el de 
Hasofer, con una variación del índice de 
confiabilidad Af = 0.3113% y del 0.1% en 
cuanto a la probabilidad de falla P,.. 


4. El siguiente resultado en aproximación es 


el que se obtiene con el método de Harr 
con AB = 3.76% y del 13.79% en cuanto a la 
probabilidad de falla P.. 


5. Los resultados obtenidos con el método de 


la transformada de Mellin y los obtenidos 
con el método de Rosenblueth son muy 
similares (difieren entre sí un Af = 0.581% 
en cuanto al índice de confiabilidad y un 
2.329% en cuanto a la probabilidad de 
falla). 


6. El método FOVE y el de primer orden (que 


se basan en las misma hipótesis) es el más 
lejano, con una variación del AP = 9.086%. 


7. Los resultado obtenidos con el método 


de Tang, por ser un análisis hecho 
en una condición límite del espacio 
multidimensional, parecen exagerados, ya 
que la variación es del AB = 65.42% superior 
al de integración directa; sin embargo, por 
ser precisamente obtenidos en la superficie 
de falla parecieran ser más realistas en 


el momento que se consideran todas las 
variables que intervienen en el análisis. 


Conclusiones y recomendaciones 


El riesgo, la confiabilidad y el análisis de 
incertidumbres pueden aplicarse en proyectos 
de ingeniería hidráulica, al estimar la proba- 
bilidad de falla en cuanto a capacidad 
hidráulica de alcantarillas y conducciones, 
desbordamiento, filtración, sismos, etcétera, en 
presas y diversos problemas más. 

Cada método tiene sus propias hipótesis 
y limitaciones, ventajas y desventajas. La 
mayoría de los métodos puede aplicarse y 
es posible obtener estimaciones confiables 
cuando los métodos se aproximan a 
comportamientos lineales, pero cuando se 
tiene una no linealidad de las variables o bien 
cuando las incertidumbres se incrementan 
significativamente, la precisión de algunos 
métodos se deteriora rápidamente. Tal es el 
caso del método del primer orden del segundo 
momento estadístico. 

Para métodos con muestras de gran tamaño 
en sus variables, el método de Monte Carlo es 
el de mayor aplicación, pero su confiabilidad 
converge cuando se tiene un gran número de 
simulaciones y no se conoce estrictamente 
el resultado final de la probabilidad de falla; 


otra limitante importante es que el número 


Cuadro 11. Resumen de resultados obtenidos con los distintos métodos. 


Metodo de cvaliáción Media Varianza Confiabilidad a Riesgo de falla 
[lo Var (W) B ». Ds 
Integración directa 0.15853 0.0059765 2.0506 0.97818 0.02182 
Monte Carlo - - - 0.97967 0.02011 
Transformada de Mellin 6.14 10.144225 1.9278 0.973197 0.0268 
FOVE 5.96 10.05172 1.8798 0.96990 0.0301 
Rosenblueth 6.22 10.022 1.939 0.97381 0.02619 
Harr 6.143 9.6612 1.9763 0.97562 0.02438 
MFOSM (primer orden) 5.96 10.05172 1.8798 0.9699 0.0301 
AFOSM (Hasofer) 5.536 7.241481 2.057 0.9802 0.01983 
AFOSM (Tang) - - 3.39211 0.99965 0.000349 
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de variables puede hacer que el problema no 
tenga una solución práctica. 

Métodos en los que se emplea el punto de 
estimación (Rosenblueth y Harr), por ejemplo, 
pueden ser muy atractivos desde el punto 
de vista computacional, en la medida en la 
que el número de variables se incrementa y 
pueden parecer muy buenos en su empleo, 
ya que ofrecen resultados parecidos a los 
obtenidos con el método de Monte Carlo y el 
de integración directa; sin embargo, en caso 
de que las incertidumbres sean importantes, 
pueden existir diferencias significativas con 
tales métodos. 

El método del primer orden del segundo 
momento estadístico es aplicable sólo en 
casos muy sencillos, en los que la función 
de comportamiento está claramente definida 
y existe una linealidad en las variables; sin 
embargo, en problemas complejos pierde 
precisión rápidamente. El método del primer 
orden del segundo momento estadístico 
es muy aplicable y puede tomar en cuenta 
incertidumbres en el caso de que el analista 
decida hacer correlaciones de las variables 
que intervienen en el problema y que la 
mayor parte de las veces se asocia con 
incertidumbres; ésta parece ser una gran 
ventaja sobre los demás métodos, ya que es 
posible involucrar variables que en muchas 
ocasiones se ignoran O desprecian por no 
poder analizarlas. 

El método de Hasofer en problemas simples 
parece ser bastante apropiado; sin embargo, 
el método de Tang resulta muy atractivo al 
analizar el estado límite de falla del problema 
estudiado. 


Notación 


Los siguientes símbolos se utilizan en este 
documento: 


C(X) = matriz de covariancia de las variables 
aleatorias X. 


D = diámetro de la conducción. 

D=matriz diagonal queinvolucra las variancias 
de las variables aleatorias. 

DPA = distribución de probabilidad acumu- 
lada. 

E(W) = esperanza del modelo W. 

F, (0) = función de distribución de probabilidad 
del margen de seguridad. 

FS = capacidad de resistencia del sistema. 
f(X) = función de distribución de probabili- 
dad de la variable aleatoria X. 

MS = capacidad de resistencia del sistema. 
M,(S) = transformada de Mellin de la función 
1%). 

N = número de variables aleatorias X. 

n = coeficiente de rugosidad de Manning. 

p,= probabilidad de falla. 

p. = probabilidad de seguridad (no falla) en el 
que la resistencia del sistema excede la carga. 
Q = gasto o capacidad de descarga. 

R = riesgo. 

S = pendiente de plantilla. 

s, = vector de sensibilidad de los coeficientes 
de la función de comportamiento W(X) 
evaluado en el punto de expansión X.. 

T. = periodo de retorno. 

W(X) = función de comportamiento o desarrollo 
de las variables aleatorias X. 

X = capacidad de resistencia del sistema. 

x'= variable estandarizada. 

Y = 
sistema. 


capacidad de demanda o carga del 


Z = variable aleatoria continua. 

Q., = derivada direccional para la variable ¡$'"“ 
en el punto x.. 

B = coeficiente o índice de confiabilidad. 

e = representa los términos de orden superior 
de la serie de Taylor desarrollada para la 
función g(x). 

U = media de las variables. 

p = coeficiente de correlación entre variables 
aleatorias. 

G = desviación estándar. 

0? = varianza. 

V = matriz de eigen vectores. 

A = 
valores. 


matriz diagonal compuesta de eigen 
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Anexo I. Transformada de Mellin 


Similar a la convolución de las transformadas 
exponencial y de Fourier, la transformada de 
Mellin se puede asociar con múltiples variables 
estocásticas independientes, en un producto 
que es simplemente igual al producto de las 
transformadas individuales de cada variable, 
asociada con su función de distribución de 
probabilidad. 

Los cuadros A.1 y A.2 muestran la trans- 
formada de Mellin de algunas distribuciones 
comúnmente usadas. 

El cuadro A.3 muestra la transformada de 
Mellin de algunas funciones de densidad de 
probabilidad comunes. 
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Propiedad E Variable aleatoria Transformada de Mellin 
Estándar Fx) XxX Ms) 
Escala fax) XxX azM (s) 
Lineal af (x) XxX aM (s) 
Translación x f(x) XxX Mi. (s+a) 
Exponenciación fx) XxX a M (s/a) 


Fuente: Park (1987). 


Cuadro A.2. Productos y cocientes de la transformada de Mellin sobre variables aleatorias estocásticas. 


Variable aleatoria FDP obtenida 
W=X fx) 
W=X fx) 

W=1/X fx) 
WE 07 ASES) 
W=X/Y FLO, 2 Ly) 

W=axoY ELO, 8, 4y) 


Ms) 
Ms) 
Mbs — b +1) 
MyQ -s) 
M,(s) M,(s) 
M,(s) M,(2 —s) 


ariM (bs —b +1) M,(cs-c+1) 


Fuente: Park (1987), a, b, c: constantes; X, Y, Z: variables aleatorias. 
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Cuadro A.3. Transformada de Mellin para algunas funciones de distribución de probabilidad (FDP) y 


de densidad f(x) comúnmente usadas. 


Probabilidad FDP, fx) Transformada de Mellin 
1 ba? 
Unif 
niforme e id 
Nao 1 ¿en 29r 5) 
ala ca 21 hb para s impar, o para s par 
Jr 
1 Lin XA x 2 1 di 
Log normal ra 2-9 exp |(s—1) yn x dei O nx 
TUXO o x 
Exponencial Pe?" p=T(s) 
B al, -px p="T(a +S -1) 
Gamma —— (pode A 
(o) P T(o) 
= a <x<m 
NA s Ss E 
Triangular 2 b(0'=m) _ alma) 
2 E ip s(s+1) (b-a)| b=m m-a 
b-a|b=m 
a x-é E x-É i O (sl largsirp(K 
j - ex ss —IBPEFTI|= +1 
Weibull 5 p (55) ejes l ) 
CCA Pla +D)a+s=) 
Beta no estándar Blo, BNO a) T(a)Tla+b+s-1) 
3(s-1 .1: : 
Pay joo 
Beta estándar ap). ol K 
a, 
donde M,(k) para la beta estándar 
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Abstract 


MARENGO, H., ARREGUÍN, El. é: ROMERO, 1. Risk assessment in hydraulic engineering 
projects: uncertainties and reliability. Water Technology and Sciences, formerly Hydraulic 
engineering in Mexico (in Spanish). Vol. I, No. 4, October-December, 2010, pp. 5-35. 


This article presents the procedure and application of several methods of risk and reliability 
analysis to a simple problem, such as the discharge capacity of a sewer. The analysis considers 
that the three variables ofroughness, diameter, and slope have normal probability distributions 
and are applied for comparing direct integration methods; the Monte Carlo Method; Mellin 
Transform, with first-order variance estimation; Rosenblueth's and Harr's Point Estimation 
Method; the First Order Second Moment Method (MPOSM); and two versions ofthe Advanced 
First Order Second Moment Method (AFOSM): Hasofer-Lind's and Tang's. The result of the 
probability of failure obtained with the direct integration method is taken as true (which is 
not possible in complex analyses, especially when uncertainties are significant) and results 
are compared. Most methods can be applied, and it is possible to obtain reliable estimations 
when the methods come close to linear behaviors, but when the variables are not linear or 
when uncertainties increase significantly, the accuracy of some methods deteriorates rapidly. 
Such is the case of the MFOSM method. For methods with samples with sizeable variables, 
the Monte Carlo method is the most commonly applied, but the reliability of the method 
converges when there is a large number ofsimulations, and the final result of the probability of 
failure is strictly unknown; another important limitation is that the number of variables may 
cause the problem not to have a practical solution. Methods where point estimation is used 
(Rosenblueth and Harr) may be very attractive from a computational perspective in as much 
as the number of variables increases and may seem to be good to use, since they offer similar 
results to those obtained with the Monte Carlo method and the direct integration method. 
However, if uncertainties are significant, there may be meaningful differences. The MFOSM 
method is applicable only in very simple cases where the behavior function is clearly defined 
and there is variable linearity; however, it rapidly looses accuracy in complex problems. The 
MFOSM method is quite applicable and can take into account uncertainties in case the analyst 
decides to make correlations of the variables that intervene in the problem and that most of 
the times are associated with uncertainties. This seems to be a great advantage over the other 
methods, since it is possible to involve variables that are many times ignored or undervalued 
because they cannot be analyzed. The Hasofer Method seems to be quite appropriate for simple 
problems; however, the Tang Method is quite attractive for analyzing the failure limit state of 
the problem studied. 


Keywords: risk assessment, uncertainty, reliability. 
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